Generalized Einstein-Maxwell field equations in the Palatini formalism 
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We derive a new set of field equations within the framework of the Palatini formalism. These 
equations are a natural generalization of the Einstein-Maxwell equations which arise by adding a 
function J-{Q), with Q = F a/3 F a p to the Palatini Lagrangian f(R, Q).The result we obtain can be 
viewed as the coupling of gravity with a nonlinear extension of the electromagnetic field. In addition, a 
new method is introduced to solve the algebraic equation associated to the Ricci tensor. 
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I. INTRODUCTION. PALATINI 
/(i?, Q)THEORIES REVISITED 

In the last few years, different versions of modified theo- 
ries of gravity have experimented wide interest in the lit- 
erature. In particular, f(R) Palatini theories have been 
regarded among the most promisin g at temps to general- 
ize Einstein's theory of gravity |^,|T^].(l3[.(l3[.p^.This 
is because of several reasons. First of all, this formalism 
provides a very elegant way to derive a cosmological con- 
stant and therefore a possibility to explain the observed 
cosmic speedup. On the other hand, it exists an intrinsic 
theoretical interest in the extensions of general relativity 
(GR),when it comes to incorporate quantum mechanics 
to a theory of gravity. 

However, not only to fully capture the phenomenology 
of GR, but also for being able to describe a richer va- 
riety of physical effects, we should go beyond f(R) the- 
ories. We can obtain this within the Palatini approach 
by just enlarging the class of Lagrangians from f(R) to 
f(R,Q), where Q = R^R^. 

The field equations for these class of Lagrangians are, for 
the metric and the connection respectively 



gians are 



JrR^u 



V a [^-g{f R g^ + 2f Q RPr)] ={ ) 



(1) 



(2) 



with f R = d R f, and f Q = d Q f. 

The general algorithm to attack these equations was 
described fn[l|], and consists in several steps. First, we 
need to find a relation between R^ v and the matter 
sources. Rewritting ([1]), using P^ — R l _ ia g aL ' we find, 



connection can thus be expressed as the Levi-Civita con- 
nection of an auxiliary metric which is related with 
9fiv by a non-conformal relation [1] 

The solution of ([3]) is only known in some particular 
cases, like the perfect fluid or the scalar field, but we point 
out here for theoretical purposes that a general solution 
exists and can be found explicitly. In order to prove this 
point, let us rewrite ([3]) in matrix notation 



VqP 2 + fn.P - -fl = k 2 f 



(4) 



wich has the structure of a quadratic matrix equation, 
AX 2 + BX + C = 0. 

We can make the following identification, 

A = 2f Q I, B = f R I, C = + fc 2 f ) 

It is straightforward to show that the quadratic matrix 
equation (4) , satisfies all the conditions required [J. This 
allow us to complete the square in the usual way to find 
the following solution for /g ^0 



P{f) 



h 



4/q' 

where we have defined 
1 



2/q 



al + pT 



«= j(/I + 4/-/q) 



P = 2fc 2 f Q 



(5) 



(6) 



(7) 



(3) 



Therefore, if the matrix al + /3T, has a square root, an 
explicit solution P(T) will always exist, but the mathe- 
matical tools required to compute the square root of this 
matrix will depend on the particular problem chosen. 



This can be seen as a matrix equation, which estab- 
lishes an algebraic relation P^ = P^{T^). Once the so- 
lution of ([3j is known, ([2]) can be written is terms of 
g^ LV and the matter, which allows to find a solution for 
the connection by means of algebraic manipulations. The 
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A. Solving for a diagonal matrix 

When a I + (3T is diagonal the solution will be auto- 
matic. We are going to see now an example of this type, 
and in the next subsection we will present a general al- 
gorithm to proceed in other, more general, situations. 
The energy momentum tensor of a perfect fluid can be 
written as 



T a p = (p + p)u a Uf3 + pg aj 3 



(8) 



where p is the pressure of the fluid and p its den- 
sity. Writing the last equation in matrix notation, after 
a bit of algebra we find the following result 



at + ftf = ( a -JP | ) (9) 
1 (a + (3p)I 3x3 



which is a diagonal matrix, and therefore computing its 
square root will be straightforward: 



aI + f3T 



yja - j5p 



(10) 



In the last expression, it was selected the positive sign 
of the square roots of the coefficients in order to be con- 
sistent with the limit Jq — > 0. These results allow us to 
write (5) for a perfect fluid as 



P{f) 



'n 7f N 
, o uji 3x3 



where 



CI 



4/q 



2^rm> - fa 

4/q 



(11) 



(12) 



(13) 



We should point out that all these results are consis- 
tent with the formalism developed in [lj , , [l^|for the 
perfect fluid, but the method described here provides a 
powerful and direct computation of the matrix P(T), a 
calculation that in some particular cases may be almost 
automatic. 



If S € TL 4x4 , then there exists a real orthogonal matrix 
U, such that 



u 1 su = w 



(14) 



where W is an upper dia gon al matrix, known in linear 
algebra as the Schur form[ll(,of the matrix S. If we can 
find a matrix L which is a square root of W (L 2 = W), 
then it can be proved that the matrix 



J = ULU 1 



(15) 



satisfies the identity J 2 = S, and therefore is a square 
root of al + f3T. In order to describe the several steps 
needed to complete the algorithm, we have to make some 
clarifications. Schur's theorem guarantees that given any 
square, real matrix S, the descomposition (14) exists. We 
first need to compute the eigenvalues and eigenstates of 
S, and then we can construct the matrix U, builded by 
eigenvectors of S. The next step requires to compute 
a square root L of the upper diagonal matrix W. To 
achieve this goal, we can make use of the relation L 2 = 
W, to write for j >i, 



(16) 



k=l 



This equation can be descomposed into two other equa- 
tions. First we have 



(17) 



and for j > i, 



3-1 



LuLij + LijLjj — Wij — y LikLjk (18) 



k=i+l 



Therefore, equation (18) provides an algorithm for 
computing the remaining blocks Lij once known the di- 
agonal blocks La. The final step will be the computation 
of the square root of 5* by means of the transformation, 
J = ULU T . 



B. A method based in the Schur descomposition 
for the general case 

Here, we will attempt to provide an algorithm to com- 
pute the square root of (5) in the general case. Let us as- 
sume that the matrix S = al + /3T is nonsingular, which 
means that its determinant does not vanish. 
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II. f(R, Q) LAGRANGIANS WITH A FUNCTION 



Inserting this last result together with the variation of 
the determinant of the metric in (23) we get 



The Palatini action that leads to the field equations 
([I]) and ([2]) for the metric and the connection is 



%,r,?Am] = -£2 J d A x^f~^gf{R 1 Q) 



+ %,</>»] (19) 



where g a p is the space-time metric, T^o is the connec- 
tion, which is indepentent of the metric, tj) m represents 
the matter fields, and finally, the scalars R = g^Rv , 
Q = R a/3 R a P. We want to consider here the possibility 
of replacing in the same action the term S\g,ip m ] by a 
function of scalars of the field strength . Taking into 
account the identity g^F 1 *" = 0, due to the fact that F^ v 
is skewsymmetric, the most natural choice is consider a 
function F(Q), with Q = F a ^F afi 

With this replacement, we find convenient to write the 
action as follows: 



S\g,T,A] 



1 



-g[f(R,Q)-k'F(Q)\ (20) 



where is a second independent connection, defined 
by means of F^ v — V ^A v — V^A^, and k' is a con- 
stant with the appropiate dimensions. In general, the field 
strength F^ u in terms of T and the potential vector A will 
be 



V„A„ - V„A M = d^A v - d v A, L - (I* - TP)A P (21) 



In the following we will assume 



T£ M which 



means that we set the torsion to zero. 
In these conditions the last term of (21) vanishes and 
therefore, V \ l A u - V V A^ = d^A v - d^A^. 

We must point out that (20) represents a particular 
case of Nonlinear Electrodynamics (NED) coupled to 
gravity.A family of f(R) Palatini Lagrangians coupled 
to NED have been studied recently in [16| The variation 
of the part of the action which contains the term f(R, Q) 
can be found elsewhere [l[ ,and therefore our task will be 
restricted to compute the variation of the second part 
which we will denote as follows 



-k' 



2k 2 



This variation gives the following result 



5Sa = 



k' 



d 4 x 



2k 2 

Since Q = g^ a g ufi F^ v F a p, it is easy to see that 
SQ = 2F ilOL F«8g» v + 2F^8F IMJ 



(22) 



(23) 



(24) 



SSo = —. 



k' 



,Q = '2k 2 / d x ^~ 9 

+ Tq {^F^SFpv + 2F fia F^8g^)] (25) 

where Tq = |f . 

The next step requires to express SF^ in terms of SA^ 
which can be done using the definition of F^, 



SF l „, = V ll (SA v )-^ v (8A li ) 



(26) 



Using the last identity in (25), we obtain after grouping 
terms: 



SS Q = -C I d l xy/=g[[--g^ + 2F Q F^ a F^)5gi a ' 



-1TqF» v {V^SA^-VuiSA^))] 



(27) 



k' 



with C = — ^-.It is convenient now to split the action 
2k 

into two parts in order to compute properly the last con- 
tribution 



M 



-2C J d 4 x^T Q F' w [V„(6A V ) - Vu(SA„)] (28) 



Using integration by parts and removing the total 
derivatives we find 



M = -2C / d 



(29) 



We have succeeded in expressing the variation of the ac- 
tion Sq in terms of variations of g^ v and A^. This is 
equivalent to the algorithm used in [l[ to express the 
variation of the action (19) which contains the f(R, Q) 
Lagrangian, purely in terms of varations of g M „ and the 
independent connection T^ v . 

It is possible to write (29) in a more compact form 
rearranging indices 



M 



2C / <FxV v 



-gT Q (F^S^5;F^)]SA x 



(30) 

With this result, the final expression for SSq in (27) 
will be: 



SSq = -C J d 4 x^=g-[( - ^g^„ + 2F Q F^F?)6g^ 
+V„ [V=£F e (F^S* - 6»F>*) ] SA X ] (31) 
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To obtain our field equations we need to compute the 
variation of the other contribution to the action in 
(20). But as we said before, this task has already been 
done elsewhere [l[ , [l] ,and the details will be omitted 
here. The variation gives: 



SS[g,F] = — 2 



WV-gifRRiu' - 



+2f Q R^RZ)+Vp[^gK XlJ Wt] (32) 

where = f R g» v + 2f Q R^. Let us define now the 
following tensor 



(33) 



Note that when T = Q = F a ^F a p it would imply that 
J-q = 1, and in this case (|33p .reduces to the usual defini- 
tion of the energy-momentum tensor of the electromag- 
netic field. It is also important to note that in general this 
tensor is not traceless: 



yielding modified dynamics. Finally, combining (32) with 
all our results, we get: 



fiiR[ii> — -^fg^v + 2fQR fia R" — —T^u (35) 



Air 



V Q W=g{fR9 M + VqR 01 )] = (36) 



VpW=gT a FW] = Q (37) 

This is a system which reduces to the usual Einstein- 
Maxwell field equations of GR when f(R, Q) = R, and 

HQ) = Q- 



I I uv 



(34) 
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